Dimer diffusion in a washboard potential 
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The transport of a dimer, consisting of two Brownian particles bounded by a harmonic potential, 
moving on a periodic substrate is investigated both numerically and analytically. The mobility 
and diffusion of the dimer center of mass present distinct properties when compared with those 
of a monomer under the same transport conditions. Both the average current and the diffusion 
coefficient are found to be complicated non-monotonic functions of the driving force. The influence 
of dimer equilibrium length, coupling strength and damping constant on the dimer transport 
properties are also examined in detail. 
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I. INTRODUCTION 

One particular example of Brownian motion on a peri- 
odic substrate is the diffusion of atoms and molecules on 
crystal surfaces This mechanism is of both concep- 
tual and technological interest Q , being relevant to het- 
erogeneous nucleation, catalysis, surface coating, thin- 
film growth, etc. Individual atoms diffusing on a surface 
can eventually meet and form dimers or trimers. For 
example, on the semiconductor Si(lOO) or Ge(lOO) sur- 
face, most of the deposited Si or Ge atoms form dimers. 
Atoms adsorbed on metal surfaces may also form closely 
packed islands that diffuse as a whole [1, 0, @. This 
raises the issue of the role of the internal degrees of free- 
dom on the transport of extended objects through micro- 
and submicro-devices. 

One of the most important problems in modern nan- 
otechnology is how to manipulate small particles in or- 
der to perform a preassigned operation. For instance, 
the mobility and diffusivity of atoms adsorbed onto crys- 
tal surfaces can be controlled by applying deterministic 
forces @, 0]- A direct manipulation method consists in 
applying a constant direct current (dc) local electric field 
by means of a scanning tunnel microscope tip Q. A se- 
lected adatom or admolecule with nonzero charge will 
then move in the direction of the electric force; neutral 
particles will be forced into a region of a stronger field 
due to induced polarization Q. This problem can be 
modeled as a Brownian motion on a tilted periodic two- 
dimensional (2D) substrate. 

In this work we study the transport of a dimer con- 
fined on a periodic substrate with a focus on the effects 
of the internal degrees of freedom on its mobility and dif- 
fusivity. For simplicity, we restrict our analysis to sub- 
strates in two or higher dimensions, which can be effec- 
tively reduced to one-dimensional (ID) systems. In the 
simple case of a dimer driven by a constant force oriented 
along a symmetry axis of a 2D substrate, one wants to 
characterize the stationary transport in the force direc- 
tion, whereas transverse diffusion is not affected by the 
bias; for a full 2D treatment, see, e.g., Ref. Of 
course, the results of the present paper apply well also 



to a variety of physical and biological systems, where the 
particle dynamics is naturally constrained to (quasi-) ID 
substrates. Examples of current interest include colloids 
or cold atoms [ijj in optical traps, superconducting 
vortices in lithographed tracks [l3| , ion-channels [Till , cell 
membranes [l5j |. artificial and natural nanopores [16| |. etc. 

This paper is organized as follows: In Sec. [TT] we in- 
troduce the model, define the units and give the details 
of our numerical simulations. Our numerical results are 
presented in Sec. IIIIl In particular, the role of the dimer 
length in the transport properties is studied in Sec. IIII Al 
the monomer like regimes (for weak and strong couplings) 
are discussed in Sec. IIII Bl finally, the influence of the 
coupling strength and of the damping constant on the 
dimer transport are analyzed in Sec. IIII CI Potential 
applications of our results to ID irreducible devices are 
sketched in Sec. IIVI 



II. MODEL 

A monomer moving on a ID periodic substrate with 
potential Uq(x) — Uq(x + L) under the influence of an 
external dc bias F and at finite temperature T can be 
described by the Langevin equation (LE), 

dUo(x 



-TjX 



dx 



+ F + €(t) . 



(1) 



Here r\ = rwy is the viscous friction coefficient, with 7 be- 
ing a damping constant and m the mass of the Brownian 
particle. The stochastic force £(i) represents the environ- 
mental fluctuations and is modeled by a Gaussian white 
noise with zero mean, (£(t)) = 0, and auto-correlation 
function (ACF), 

= 2r,k B T8(t-lf). (2) 

For a symmetric dimer the corresponding LE's have 
the form 
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where £i(t), i = 1,2, are two independent zero-mean 
stochastic processes with ACF, 



(Zi(t)Z j (t'))=2r ] kBT5 ij 5(t-t'). 



(4) 



Note that the inter-particle interaction is incorporated in 
the potential function, 

K 

U(x 1 ,x 2 ) = U Q {x x ) + U (x 2 ) + — (x 2 -X\- a ) 2 ■ (5) 

That is, we assume the interaction between the two dimer 
particles to be harmonic with coupling constant K and 
equilibrium distance clq. The simplest choice for the pe- 
riodic substrate potential is [l7l ]. 



Uq(x) = Aq cos(kx) . 



(6) 



with k = 2tt I L. 

The LE's ([I]) and ([3]) can be conveniently rescaled. By 
introducing suitable space, energy, and time units, 



(7) 



X=l/k, e = Aq , t = \J \ 2 m/e , 
we define the dimcnsionlcss quantities: 

x . a ~ k B T ~ A 

X =T, a Q = —, T= , F = -F, 

A A e e 

K=-K, t = -, 7 = 7T, l(i) = -m- (8) 
e r e 

No particle can be trapped by the potential (J6J) under 
any circumstances for tilting larger than the critical value 
F cr = 1 (in rescaled units). In the following we drop the 
tilde altogether. 

After rescaling, the LE (p} for a monomer moving in 
the potential ^ reads, 



= -jx + sinx + F + £(t) 



(9) 



where the ACF of the rescaled noise is £,(t')) — 
2jT6(t - t'). Analogously, the coupled LE's © for a 
symmetric harmonic dimer in the same substrate poten- 
tial become [see Eq. ([5])], 

xi = -7.T1 + sinxi + F + K(x 2 — x\ — a ) + i\{t) , 
x 2 = -7.T2 + sina;2 + F — K(x 2 — x\ — a ) + ^(t) , 

(10) 

with (*0) = 2 1 T5 lj 5{t - t'). 

The dimensionless LE © for a monomer and (fTO]) for a 
dimer have been integrated numerically through a stan- 
dard Milstein algorithm [l^|. Individual stochastic tra- 
jectories were simulated for different time lengths i max 
and time steps At, so as to ensure appropriate numeri- 
cal accuracy. Average quantities have been obtained as 
ensemble averages over 10 4 trajectories; transient effects 
have been estimated and subtracted. 



III. RESULTS: MOBILITY AND DIFFUSION 

When considering a pair of interacting Brownian par- 
ticles, it is natural to study the motion of their center of 
mass, 



X = i(lEi + X 2 ) . 



(11) 



The quantities that best characterize the stationary 
dimer flow are: (a) the net velocity, 



v = lim 



(*(*)> 



(12) 



or, equivalently, the related mobility, fi — v/F; (b) the 
diffusion coefficient, 



(13) 



where (6X 2 ) is the mean square displacement of the cen- 
ter of mass, i.e., 



(SX 2 ) = (X 2 ) — 



(14) 



- 4(8x1) + \(8x 2 2 ) + ~{{x lX2 ) - (xi)(x 2 )). 

For the following discussion we also introduce the rel- 
ative coordinate F, 



Y = x 2 - xi , 



(15) 



representing the dimer size. The quantity Y can in prin- 
ciple also become also negative. However, this happens 
only when the dimer oscillations around the equilibrium 
position become very large. In the range of parameters 
adopted in the present paper, we have verified that Y 
remains positive even for small values of the elastic con- 
stant K, where one recovers the monomer limit. In fact, 
the distance Y can become negative if both monomers 
fall into the same valley. In our simulations, the dimer 
length (at rest) varies in the range do € [L,2L]. Thus, 
the monomers start out in different potential valleys and 
are observed to stay so for all times (i.e., configurations 
with Y < do not occur). 

The LE's © can be rewritten as a LE for the center of 
mass coordinate X and one for the dimer length Y, that 



X = - 7 A + cos(r/2)sinA + F + Q(t)/V2, (16) 
Y = -7y + 2cosAsin(y/2)-2^(y-a ) + \/2g(t). 

(17) 

Note that the two noises Q(t) = [^(t) + &(t)]/V2 and 
li't) — ~ ^iW]/v^2 are uncorrelated and have the 

same statistics as £,i }2 (t), namely, (q(t)) = (Q(t)) = 
and 



(«(<MO) = (Q(t)Q(t'))=2 7 TS(t-t') 



(18) 
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Figure 1: (Color online) Mobility (a) and diffusion coefficient 
(b) versus the dimer length ao for different values of the tilting 
force F. Simulation parameters: coupling constant K — 1.5, 
temperature T — 0.1 and 7 = f . is the free dimer diffu- 
sion coefficient Doo = Dq(T/2); see text. 



In the absence of a substrate potential the mobility of 
both a monomer and a dimer is /io = 1/7. Corre- 
spondingly, the free diffusion coefficient for a monomer, 
Dq(T) = T/7, is twice as large as that for a dimer, 
D {T/2). 



A. The role of the dimer length 

At variance with a monomer, a dimer has two degrees 
of freedom. This affects its diffusion dynamics [19| to 
the point that its diffusion coefficient D can develop a 
non-monotonic dependence on the dimer parameters. For 
instance, dimer transport strongly depends on the ratio 
between the period L of the substrate and the natural 
length a of the dimer [U El HI . 

In the absence of an external force, F = 0, at low 
temperature the diffusion coefficient of a rigid dimer de- 
creases monotonically on raising the dimer length ao 
from L/2 to L. This can be well understood from 
Eq. p^|) . In the limit K = 00 the dimer length is just 
Y = ao and the force cos (Y/2) sinX acting on X(t), 
Eq. (fTT5|) , corresponds to a periodic potential with ampli- 
tude I cos (Y/2)\. For ao = L/2 — it this quantity is zero 
and the dimer center of mass undergoes free diffusion. 




Figure 2: (Color online) Dimer configurations corresponding 
to zero pinning force and maximum diffusion coefficient; see 
also Fig. [T] and text. 



For ao = L = 2tt the periodic potential amplitude is 
maximum, | cos(Y/2)| = 1; diffusion in a periodic poten- 
tial is known to be suppressed compared to free diffusion 
[23l [23|. Therefore, the maxima and minima of D versus 
ao coincide with the minima and the maxima of the mod- 
ulating factor |cos(ao/2)|, respectively. This conclusion 
applies also to the case of finite elastic constants as long 
as (Y(t)) w ao, that is for rigid dimers, K 3> 1, at low 
temperatures, I < 1. (For the opposite limit of weak 
dimers, K < 1, see Sec. IIIIBI ) 

In the presence of a sub-threshold external force, F < 
F CI , the diffusion coefficient D is a nonmonotonic func- 
tion of the dimer length ao, as shown in Fig. [ljb). The 
numerical results in Fig. [T] have been obtained by sim- 
ulating a relatively rigid, K = 1.5, and moderately 
damped, 7 = 1, dimer. In the case of a strong to mod- 
erately damped monomer in a washboard potential, the 
curves D(F , T) are known to develop a peak around F cx , 
where the barrier height of the tilted periodic potential 
Uo(x) — Fx vanishes [25(. Analogously, in the case of 
a dimer, D attains a maximum for dimer lengths such 
that the effective pinning force also vanishes, i.e., for ao 
equal to the distances between maxima and minima of 
the washboard potential, see Fig. [5] In the case of a 
driven rigid dimer with F < F cr = 1, this takes place 
for equilibrium lengths — (L/2)[l± (2/ir) arcsin(F)]. 
Note that are given mod(L) and 0,^+0,^= L. 

Figure [Ua) demonstrates that the mobility is small- 
est for commensurate dimers with ao = L and largest 
for ao = L/2 (see also Ref. 22]). The smaller the ap- 
plied constant force, the smaller is the ao range around 
ao = L/2, where the mobility of the dimer is significantly 
different from zero. For large enough tilting the dimer is 
considerably mobile, no matter what is the value of ao- 
For F —> 00 the mobility fj, — > [1q and the effective dif- 
fusion coefficient D — > Doo = Dq(T/2). We remark that 
the ao dependencies of fj, and D shown in Fig. [T] are given 
mod(L) [22j |. In fact, the system dynamics, as given by 
Eqs. (|10p . is invariant under the change ao — > oq + L and 
%l —> %1 — L (or X2 — > X2 + L). 
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Figure 3: (Color online) Diffusion coefficient D versus the 
tilting force F for a dimer length ao/L = 1 and different 
coupling constants K; T = 0.1 and 7 = 1. The corresponding 
curves for monomers of temperature T (solid, black) and T/2 
(solid, red) are drawn for comparison (see text). 



B. Monomer like regimes 

In the case of a monomer the mobility and the dif- 
fusion in a tilted periodic potential are in general well 
understood. In the low-temperature regime, T« 1, the 
particle mobility is close to zero for sub-threshold tilt- 
ing (locked state). Around a depinning threshold Fd, 
the mobility grows sharply and in the large force limit 
reaches the free particle limit ^o (running state). If the 
temperature is increased, the transition from the locked 
to the running state is smoother. In the overdamped 
regime, convenient fully analytical expressions are avail- 
able for both the mobility fJ,(F, T) (the Stratonovich 
formula [26|]) and the diffusion coefficient D(F,T) (the 
Cox formula [13] )• For small biases and low temper- 
ature, the diffusion coefficient is suppressed compared 
to the free diffusion Dq(T); in linear response theory 
D(F,T) ~ n(F,T)T [17|] . Depinning occurs around the 
critical tilt, i.e., F d ~ F cr , as signaled by D(F,T) over- 
shooting Do [251 ]; the lower the temperature, the more 
prominent is the growth of the depinning diffusion peak. 
In the large force limit, the free diffusion regime Dq(T) 
is eventually recovered. 

In the underdamped limit, 7 <C s/ F CI , the mobility 
and the diffusion coefficient display a similar behavior 
with one significant difference: the depinning threshold 
Fd is a monotonic function of the damping constant with 

lim F d ~ 3.36 7 v/^ (19) 

and F d ~ T cr for 7 > ^W, |H S3. 

In the case of a dimer, the general behavior recalls that 
of a monomer, namely, both the transition of the rescaled 
mobility from to [1q and the corresponding enhance- 



ment of the diffusion coefficient above its free diffusion 
value still occur as the tilting force is increased past the 
depinning threshold. The monomer dynamics is a use- 
ful benchmark to check the accuracy of our simulations 
for the dimer diffusion. Indeed, in the limit K — > 0, 
Eq. US]) boils down to (5X 2 ) = (5x\)/2, with xi obey- 
ing the monomer LE ((5) with temperature T. It follows 
that for a weak dimer, K <C 1, the ratio D/D^ is closely 
reproduced by the analytical curve D(F,T) / Dq(T) ob- 
tained from the monomer LE @. This argument ap- 
plies to both commensurate, Fig. [3J and incommensurate 
dimers, Fig. QJb). 

Rigid dimers also behave like monomers. In the limit 
K — ► 00, the solution of Eq. |T7) is Y(t) = Y = a 
and Eq. (fT6|) is then equivalent to the monomer LE (J9j) 
with temperature T/2 and substrate amplitude (critical 
tilt) cos(ao/2). Accordingly, for commensurate dimers 
with ao equal to an integer multiple of the substrate 
constant L, the ratio D/Doq is reproduced by the curve 
D(F,T/2)/D (T/2) obtained for a monomer on a tilted 
cosine potential with amplitude | cos(ao/2)| = 1 and tem- 
perature T/2 (see Fig. 13). 

Note that for large values of damping the monomer 
curve can also be computed analytically through the Cox 
formula [27| • The data in Fig. [3J confirm that for increas- 
ingly large K the depinning threshold approaches T cr = 1 
from below, as the effective critical tilt (| cos(V>/2)|) tends 
to unity. Not surprisingly, for the commensurate dimer 
of Fig. [3J the mobility curve coincides with the monomer 
mobility ^{F, T), in the weak coupling limit, and with the 
monomer mobility at half the temperature T, /i(T, T/2), 
in the strong coupling limit; both limiting curves are 
closely approximated by the Stratonovich formula (not 
shown) . 

For K — > 00 incommensurate dimers behave like 
monomers moving on a tilted cosine potential with ampli- 
tude I cos(ao/2)| < 1 and temperature T/2 (see also Fig.[S] 
for a finite coupling). When ao is equal to a half- integer 
multiple of the substrate constant L, the amplitude of 
the effective substrate acting on the dimer coordinate X 
vanishes, |cos(ao/2)| = 0, and the dimer diffusion be- 
comes insensitive to the substrate, with mobility /iq and 
diffusion coefficient Do (T/2). 

Figures [4] and [5] indicate that for a finite K the dimers 
exhibit a much more complicated behavior, which will be 
discussed in the forthcoming section. 



C. The dependence on the coupling strength 

The problem of a dimer diffusing in a washboard poten- 
tial has been studied in fact in many papers, but due to 
the large parameter space, important effects went unno- 
ticed. In Ref. [29j], it was found that for a commensurate 
dimer, D(F, T) had two maxima as a function of the tilt- 
ing force F, whereas an incommensurate dimer behaved 
more like a monomer, with D showing only one peak. 
However, as shown in Fig. (Hb), one can observe two F 
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Figure 4: (Color online) Mobility (a) and diffusion coeffi- 
cient (b) versus the tilting force F for an equilibrium distance 
ao/L — 1.5 and for different values of the coupling constant 
K; T = 0.1 and 7 = 1. In both panels the results are com- 
pared with the corresponding monomer curves (see text). 



Figure 5: (Color online) Mobility (a) and diffusion coefficient 
(b) versus the tilting force F for a coupling constant K = 1.5 
and for different values of the equilibrium distance ao. T and 
7 have the same values as in Figs. [4] and [3] 



maxima also in the diffusion coefficient of a noncommen- 
surate dimer; correspondingly, the mobility curve /z ver- 
sus F develops the nonmonotonic behavior displayed in 
Fig. BJa). More remarkably, for the same temperature 
and damping constant of Fig. 2J commensurate dimers 
presented a single peaked diffusion coefficient and mono- 
tonic mobility as functions of the tilt (see Figs. [3] and [5]). 
However, for different simulation parameters (like those 
in Ref. [2!|) two-peaked D curves were detected for com- 
mensurate dimers, as well. Thus, a doubly peaked dif- 
fusion coefficient is no signature of dimer-substrate corn- 
mensuration: the coupling constant (Fig. 2J, damping 
constant (Fig. [5]) , and temperature also play a significant 
role [see Eqs. 1221 and [25]. 

To investigate the origin of the two competing diffusion 
mechanisms shown in Fig. 2J we address in detail the case 
of a dimer with length ao equal to a half-integer multiple 
of the substrate constant L. For a finite coupling strength 
K, on setting Y(t) = a + tp{t), the coupled LE's (fl"6|) 
and p7p read 

X = -~fX -sm(i[>/2)sinX + F + Q(t)/y/2, (20) 
t/j = -~fip + 2 cos (ip/2) cos X -2Kip + V2q(t). (21) 

If the dimer is sufficiently rigid and the tilting force F 



weak, then ip(t) is small and mostly controlled by ther- 
mal noise. From Eq. (|2~lj) , on neglecting the substrate 
force with respect to the dimer coupling, energy equipar- 
tition yields (■0 2 (t)) = T/K. Moreover, the force term 
sin (i/ , /2) sin X in Eq. (f2"0"| can be treated as resulting 
from a randomly flashing cosine potential with ampli- 
tude 2(| sin [^(t)/2]|) « \tp\. This can be regarded as 
an instance of the "parametric resonance" approach pur- 
sued by the authors of Ref. 20} in the limit T = 0. On 
assuming a Gaussian distribution for tp, a corresponding 
7-independent effective critical tilt can thus be estimated, 
namely, 

Fr a [(2/Tr)(4, 2 (t))] 1/2 = V^M<- (22) 

As pointed out in Sec. IIII B[ for large to intermediate 
values of damping, the critical tilt coincides with the 
effective dimer depinning threshold Fd- For K > 0.2, 
Eq. (|2"2"|) locates rather accurately the first F peak of the 
simulated diffusion coefficient reported in Fig. SJb) . 

For F > F\ both the dimer mobility and the diffu- 
sion coefficient tend towards their free particle values, 
unless an internal resonance sets in. Indeed, driven by 
a strong force F, the dimer center of mass acquires an 
almost constant speed F/j. On inserting X(t) ~ Ft/j 
into its right hand side, Eq. (f2"Tj) becomes the LE of a 
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Figure 6: (Color online) Mobility (a) and diffusion coefficient 
(b) versus the tilting force F for ao/L — 1.5 and different 
values of damping constant 7; T — 0.1 and K — 0.3. Two- 
peaked diffusion curves are clearly distinguishable for 7 < 1 
only. 



(see Sec. MI Bp . Equivalently, incommensurate dimers 
with 7 > 2\[K must be regarded as overdamped as far as 
their internal coordinate Y is concerned; therefore, their 
diffusion coefficients are characterized by one maximum 
located around the 7-independent depinning threshold 
Fd in Eq. (|22j): see Fig.[6jb). When 7 decreases, both dif- 
fusion peaks shift towards smaller values of F. The expla- 
nation is very simple: The resonance threshold F2 tends 
almost linearly to 1/2; in the underdamped regime, the 
depinning threshold Fd is proportional to 7 as it obeys 
law (fll?)) with F cr given by the effective critical tilt F\ 
of Eq. (|2"2"|) . This estimate for Fd in the underdamped 
limit is consistent with the anticipated locked-to-running 
transition thresholds exhibited by the mobility curves of 
Fig. (Ha) with 7 < 0.3. 

Going back to the dynamics of the damped incom- 
mensurate dimer of Fig. dj we remark that on increasing 
K the resonance diffusion peaks, in addition to shifting 
to higher F (directly proportional to yK), flatten out 
on top of the plateau D = Dq(T/2); as the depinning 
peaks move to lower F (inversely proportional to \/K), 
for K — > 00 the diffusion coefficient eventually tends to 
Dq(T/2), as anticipated in the previous Sections. 

The argument presented here can be easily generalized 
to the case of commensurate dimers, or to any equilib- 
rium length; the ensuing properties of commensurate ver- 
sus noncommensurate dimers and the different monomer 
limits of the dimer dynamics have been anticipated, re- 
spectively, in Sees. IIII Al and IIIIBI 



IV. CONCLUSION 



Brownian oscillator subjected to a harmonic force with 
angular frequency £1 = F/j. Accordingly, the internal 
degree of freedom of the dimer, represented by the co- 
ordinate Y, resonates for FH approaching \J2K — 7 2 /2 
(parametric resonance [2fj|. r29l. [30 . Hl|), thus leading to 
a threshold- like enhancement of the dimer diffusion [251 ] . 
Our argument can be refined further by noticing that at 
resonance the processes X(t) and ip(t) synchronize their 
phases, so that the substrate force in Eq. (f20)) does not 
average out any more. In the presence of synchroniza- 
tion, (sin (ip/2) sinX) ~ 1/2, which amounts to replac- 
ing F with F — 1/2. In conclusion, for relatively large 
damping constants, namely 1 < 7 < 2yK, a resonance 
diffusion F peak is expected for 



F2^\+1]J2K-^. (23) 

in reasonable agreement with the simulation results of 
Fig. IUb) for 7 = 1. Correspondingly, the mobility curves 
describe a two step transition from the locked to the run- 
ning state. 

For weak dimers, K < (j/2) 2 , the two peaks of the 
diffusion coefficient tend to merge, as shown in Fig.QJb), 
and in the limit K — > a monomer dynamics is recovered 



In this paper we have studied a system consisting of 
two harmonically interacting Brownian particles diffusing 
in a ID washboard potential. We found that the aver- 
age current and the diffusion coefficient of such a dimer 
exhibit a complicated non-monotonic behavior as a func- 
tion of the driving force and the ratio of the dimer length 
to substrate constant. In the limits of the weak [K — ► 0) 
and strong (K — * 00) coupling constant the expected 
monomer dynamics was recovered. Moreover, we stud- 
ied in detail the dimer transport for different coupling 
strengths and damping constants. We concluded that 
the appearance of the second resonant peak of the diffu- 
sion coefficient versus the driving force is not related to 
the dimer length-to-substrate constant ratio, but rather 
to the damping-to-coupling constant ratio; the diffusion 
coefficient D(F) possesses two peaks only for relatively 
low damping values. 

Finally, we recall that a simple ID model is not al- 
ways a viable tool to analyze transport in two or higher 
dimensions: such a modeling makes sense for highly sym- 
metric substrates, only. There exist irreducible 2D and 
3D devices where particles are driven on an asymmetric 
potential landscape by an ac or dc driving force perpen- 
dicularly to the symmetry axis of the potential. Such a 
geometry has recently attracted broad interest [32| in the 
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context of separation of macromolecules, DNA, or even 
cells, because it is capable of inducing a transverse drift 
as a function of the drive and of the particle geometry: 
as a consequence different objects can be separated de- 
pending on their center of mass diffusion coefficient [33| . 
While the motivations of the present study apply to this 
class of devices, too, it is clear that their characterization 
must take into account the dimensionality of the system 
at hand. Dimensional reduction is limited by the spatial 
symmetry of the substrate and the particles. This is the 



subject of ongoing investigation. 
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